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Abstract 



A novel method for calculation of the motion and radiation reaction 
for the two-body problem (body plus particle, the small parameter 
m/M being the ratio of the masses) is presented. In the background 
curvature given by the Schwarzschild geometry rippled by gravita- 
tional waves, the geodesic equations insure the presence of radia- 
tion reaction also for high velocities and strong field. The method 
is generally applicable to any orbit, but radial fall is of interest due 
to the non-adiabatic regime (equality of radiation reaction and fall 
time scales), in which the particle locally and immediately reacts to 
the emitted radiation. The energy balance hypothesis is only used 
(emitted radiation equal to the variation in the kinetic energy) for 
determination of the 4-velocity via the Lagrangian and normaliza- 
tion of divergencies. The solution in time domain of the Regge- 
Wheeler-Zerilli-Moncrief radial wave equation determines the met- 
ric tensor expressing the polar perturbations, in terms of which the 
geodesic equations are written and shown herein. 
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A novel method for calculation of the motion and radiation reaction for the two- 
body problem (body plus particle, the small parameter m/M being the ratio of 
the masses) is presented. In the background curvature given by the Schwarzschild 
geometry rippled by gravitational waves, the geodesic equations insure the presence 
of radiation reaction also for high velocities and strong field. The method is gen- 
erally applicable to any orbit, but radial fall is of interest due to the non-adiabatic 
regime (equality of radiation reaction and fall time scales), in which the particle 
locally and immediately reacts to the emitted radiation. The energy balance hy- 
pothesis is only used (emitted radiation equal to the variation in the kinetic energy) 
for determination of the 4- velocity via the Lagrangian and normalization of diver- 
gencies. The solution in time domain of the Regge-Wheeler-Zcrilli-Moncrief radial 
wave equation determines the metric tensor expressing the polar perturbations, in 
terms of which the geodesic equations are written and shown herein. 

1 Introduction 

Regge and WheleeE.0 proved the vacuum stability of the Schwarzschild black 
hole, while ZerilliBH studied the emitted radiation via the radial wave equation 
for polar perturbations which source is a freely falling test mass m towards the 
black hole of large mass M. Moncriefcl showed the gauge invariant significance 
of the wave equations. Past work was concerned on the Fourier analysis of 
the emitted radiation (for a review see Ruffini til) and not on the motion of 
the particle. Radiation reaction Jor particles around black holes has been 
tackled by Cutler and co-workers u whose analysis (via Teukolsky formula) is 
based on the adiabatic approximation, the postulate of energy-balance, and 
is limited to 6 M; for_t;hc Dirac-Galts'ov method 13 antitrasformation of all 
frequencies is required u, while the axiomatic approach □ is conceived within 
the one-body problem. The aim of this work is the identification of radiation 
reaction, without the assumption of adiabacity, with minimal use of the energy 
balance postulate, and identification of the trajectory up to the horizon. Future 
developments may include a general method for motion of small objects in any 
orbit; the solution of 2nd-order equation, which energy-momentum tensor is 
based on the geodesic equations calculated herein, and a post-Schwarzschildian 
formalism. 
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2 Approach 

Past work was confined to a particle falling in an unperturbed Schwarzschild 
metric, but radiation reaction requires a geodesic path through the emission 
of radiation, in a perturbed Schwarzschild metric. To this end it appears 
necessary: a numeric or analytic solution of Regge- Wheeler- Zerilli-Moncrief 
(RWZM) equation in time domain; determination of the perturbation compo- 
nents and of lst-order Lagrangian; finally identification of radiation reaction 
by subtraction of the Oth-order terms in the geodesic equations. The RWZM 
equation for polar perturbations is: 



d 2 <f t (r,t) d 2 *;(r,t) 



dr* 2 



dt 2 



Vi(r)*i(r,t) = Si(r,t) 



(1) 



where r* — r + 2Mln ^t^J ~ ^) * s ^ e tortoise coordinate and the potential 
Vi(r) is: 



VM - (\ 2M\ 2A 2 (A + l)r 3 + 6X 2 Mr 2 + 18XM 2 r + 18M 
l{j > ~ V ~) r s (\r + 3M) 2 



(2) 



Further A = |(£ — + 2) and Si(r, t) is the 2'-pole source component and, 
for a radially falling particle, is: 



('-¥) 



(A+l)(Ar +3M) 



[r[ 1-^) S'[r-r (t)] + 



6M 



Ar + 3M 



M 



-A-l 



<5[r-r (i)] 
(3) 



where ft = iM-J (21 + l)n and ro(t) is the inverse of: 
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(4) 



Eq. (|T|) can be rewritten in terms of t,r* and with constant coefficients of the 
2-nd derivatives but solely via an approximate inverse function r(r*) and thus 
resulting into an approximate p.d.e. (further, the solution is most interesting 
at r* = oo). In the (i, r) domain instead, eq.(Q) becomes: 



1 d 2 *; 1 - A(r) dfi d 2 * ; 



where A(r 



A 2 (r) dr 2 rA 2 (r) dr 
dr* r 



dt 2 



Vj(r)tfi = S,(r,t) 



(5) 



dr r - 2M ' 



The polar perturbations are determined by: 
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where H ,H 1 ,H 2 ,h ,h 1 ,K,G are functions of t,r. The Regge- Wheeler 
gauge specifies G = ho = hi = 0. There are relations between * and derivatives 
and the perturbation components. Analytically the initial value problem is 
well defined and 4> are zero, i.e. particle at rest at infinity). Numerically, 
Lousto and Pricell3 have determined intermediate conditions, solving eq. ([[]) 
with a finite difference scheme and direct integration of the source term (thus 
an interpolation should be performed for finding * and 4/). Alternatively for 
an analytic solution of eq.(^), the asymptotic solutions for r — > oo and r — > 
would be instrumental for qualitative methods as differential inequalities or 
Lagrange identities techniqueslid. 



3 The equations of motion 

The Schwarzschild metric up to lst-order perturbations in the Regge- Wheeler 
gauge is: 



ds 2 = dsl+(l~^^H Ydt 2 +(l-^y i H2Ydr 2 +r 2 KYd9 2 +r 2 sm 2 eKYd<j} 2 +2H 1 Ydtdr 

(7) 

The geodesic equations encompass three perturbation schemes: i) a perturba- 
tive field g"" = ri% chtuarzschild v + ft^_ order perturbaUon ; ii) the particle trajectory 
r P (t) = r (t) + n(t) = r [l + e r (i)] where r (t) is the trajectory crossing the 
Schwarzschild metric and radiating out but the motion being unperturbed and 
ri(t) is the correction due to radiation reaction, via the small parameter e r (t); 
iii) the field at r p (t) is a McLaurin series at r ,t , and t p = to + ti = t [l + e t (t)\: 
g^(r p ,t p ) ~= g^(ro,t ) + e r r ^]( ro , to ) +e*to^|(r ,t ) 

There are two unknown variables e r (i) and et(t). For the case of radial 
fall d# and g?</> vanish, and using the time independence of the Schwarzschild 
metric, and dropping the (ro,io) notation one equation is: 
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The quantity d 2 r/dt 2 is deduced as ratio from eq. (||), and from the same 
eq. in where (i = t. Eq. (||) at lowest order (h = h' = h" = n" = 0) expresses the 
known motion of a particle in a unperturbed metric. Otherwise, the particle 
crosses the perturbed metric and thus its motion is influenced by the emitted 
radiation. The velocities r and t are derived from the Lagrangian. The field 
contribution at lst-order is: 
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£i =(l-—jH Yi 2 +(l-—J H 2 Yr 2 + r 2 KY6 2 + r 2 sin 2 9KY(j> 2 + 2H 1 Ytr (9) 

where the dot indicates s derivatives. Differentiating Ci (omitting the index 
I) and eliminating angular dependence, the canonical quantities are: 
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(13) 



The Hamiltonian is equal to the Lagrangian due to the absence of the 
potential. The Lagrangian is time independent for a conservative system: the 
power of the emitted radiation P gw is added to ( [To| ) pt,totai — pt + P gw = 0. The 
integration in time leads to a constant of energy dimensions from which i is 
derived: 



f ( 2M\dH „- 2 / 2M\~ 1 dH 2 ^. 2 „affi„... f „ , „ , s 

PUtotal = J (l-—J ~Qf Yt +{}-—) -Of^ +2 ^f Ytrdt + J P ^ dt 

The Lagrangian is unitary for timelike geodesies and substituting ( [l4|) in 
(^) for 9 = 4> — 0, i(r) and r are calculated. The initial conditions at infinity, 
r = and E gw — 0, fix the constant E. The zeroth order terms must be added 
to the lst-order ones. 



3. 1 Singularities 

The energy going into the black hole is E in ~ 0.3927mc 2 while the energy 
radiated away is E out ~ 0.0104/i 2 c 2 /M where the reduced mass ^ ~ m. The 
relation between E and ip or: 

T M " = 3^h;a->ST (15) 
should bring to the normalization of divergent expressions. Alternatively a 
finite size mass, e.g. dust, could be considered. 
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4 Conclusions 



The approach for identification of radiation reaction of masses falling into 
Schwarzschild black holes via lst-order polar perturbations has been shown 
and the equations of motion in symbolic form have been found. Radiation 
reaction is a fundamental concept in bodies motion theory, but also has relevant 
implications on detector's templates since the capture of stars by black holes 
is a source of gravitational waves. 
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